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Abstract 

The balance function has been proposed to 'clock' hadronization by measuring the charge- 
dependent correlation length on rapidity. According to the proposal a narrowed balance function 
would imply reduced hadron diffusion time and therefore delayed hadronization due to a long- 
lived prehadronic state, a quark-gluon plasma or QGP. I examine this chain of argument in the 
context of hadronization and rescattering. I then consider the algebraic structure of the balance 
function and its systematic dependence on the charge correlation length on rapidity. I conclude 
that the width of the balance function cannot determine a time interval from hadronization to 
kinetic decoupling, nor can the width determine the existence or temporal extent of a QGP. 

1 Introduction 

Balance function measurements were motivated by predictions [1, 2] of length-scale systematics for 
charge-dependent rapidity correlations in heavy ion collisions, according to which changes in correla- 
tion length are caused by thermal and diffusion effects marking the elapsed time from hadronization. 
Measurement of a charge-dependent correlation length on rapidity with the balance function would 
thus reveal the timing of hadronization, and detection of a reduced width of the balance function 
would indicate a long-lived prehadronic or QGP state [1]. It was suggested to vary collision central- 
ity and energy to look for such narrowing, a sudden drop in the balance-function width signaling 
QGP formation [2]. 

According to [2] both the balance function and net-charge fluctuations are sensitive to charge- 
dependent longitudinal momentum correlations. Production of a QGP of significant space-time 
extent would imply substantial production of new charge late in the collision, reducing the time 
hadron charge pairs have to diffuse apart by rescattering. A reduced balance function width or 
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'tight correlation' would thus provide evidence of a QGP. By a related argument the presence of a 
QGP would be signaled by unusually small net-charge fluctuations [3]. The theoretical hypothesis 
and balance function width interpretation depend on two aspects: 1) the relation of a long-lived 
QGP to the charge-dependent correlation length on rapidity and 2) the relation of the balance 
function width to that correlation length. The interpretation of the rms width of the balance 
function relative to the charge-dependent correlation length and the structure of the theoretical 
motivating argument are the subjects of this paper. 

I first examine the physical arguments motivating the balance function measurement, then 
review some analyses of elementary collisions which served as prototypes for the balance function, 
summarizing the methods and the understanding of hadron charge correlations which resulted. I 
restate the balance function definition in a more detailed algebraic form and thereby relate it to 
the autocorrelation distribution. I adopt a simple model of the charge-dependent autocorrelation 
as observed in heavy ion collisions and use this to investigate balance function width systematics 
relative to the charge-dependent correlation amplitude and length. 

2 Physics of the Balance Function 

According to the theoretical motivation of the balance function in [1, 2], the appearance of quark 
pairs as hadrons signals the end of the QGP phase and the start of hadron diffusion (rescattering) . 
To say that hadronization is 'late' implies that a QGP existed for an extended prior time interval. 
Diffusion alters hadron charge correlations: separation of correlated charge pairs in space increases 
monotonically with diffusion time. Schematically, the argument in [1] is: 1) Hadron charge pairs 
are created at single points in space-time; 2) thermal motion, diffusion and axial Hubble expansion 
map these points to a peaked distribution on relative rapidity, the peak width (correlation length) 
increasing with diffusion time; 3) the balance-function width measures this correlation length; thus, 
4) if the balance-function width is reduced the time from hadronization to decoupling is less and 
the time before hadronization is greater, indicating an extended QGP lifetime. 

2.1 Details of theoretical argument 

The motivating language for the balance function implies a fixed time interval between first nu- 
clear contact and kinetic decoupling, divided into two parts by hadronization. The existence and 
extent of the prior QGP interval are the nominal objects of the balance function analysis. During 
hadronization quark and/or hadron charge pairs are said to appear at nearly the same instant, and 
hadron charge pairs are created at zero space-time separation [1], which then defines the initial 
space-time autocorrelation of hadron charge pairs in this picture. 

The autocorrelation on space difference variable za maps to the autocorrelation on rapidity 
difference i/a via the phase-space correlation of axial Hubble flow or Bjorken expansion. The width 
of the corresponding structure on ?/a is determined by several factors, primarily thermal velocities 
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at kinetic decoupling and diffusion (rescattering) time from hadronization to kinetic decoupling. 
Cooling by isentropic expansion and diffusion time are both indicated by changes in the charge- 
dependent correlation length on rapidity according to the argument. 

The balance function, based on obscrvables formulated in [4], was designed to quantify these 
rapidity correlations. It is said to describe the conditional probability that a charge in one rapidity 
bin is accompanied by a charge of opposite sign in another bin. If the balance function width 
decreases (due to reduced diffusion time) the time interval prior to hadronization must have in- 
creased. 'Late' or 'delayed' hadronization inferred from the balance-function width thus implies a 
long-lived QGP. The authors of [1] conclude "...it seems clear that the canonical picture of a heavy 
ion reaction, quark-gluon plasma formation followed by late-stage hadronization, should have a 
clear signature in the balance functions." 

2.2 Critique of theoretical argument 

• Fixed total interval: The theoretical argument relies on implicit duration complementarity 
between the hadronic time interval and the prehadronic or QGP interval. Reduction of the 
former implies increase of the latter. The fixed total interval is conveyed by language like 'late- 
stage' or 'delayed' hadronization, which suggests an extended QGP lifetime. The implication 
is not justified. 

• Point-like pair creation: The assumed point-like space-time autocorrelation for oppositely- 
charged hadrons at formation is inconsistent with the fragmentation picture contained in the 
Lund model [5]. The actual hadron autocorrelation at formation should have a non-zero 
width representing QCD string physics (in the ID case), which hadronic diffusion, thermal 
velocities and other contributions then modify in the projection to rapidity difference. Rapid- 
ity correlation lengths in elementary collisions, best reflecting this base width, should provide 
a reference for A-A collisions. 

• Conditional probabilities: Analysis of measure conservation in elementary collisions [4] makes 
use of conventional conditional probabilities and generates interpretable results. The 'condi- 
tional probability' invoked in [1] is not a probability. The binning geometry is not that of a 
conditional probability corresponding to Bayes' theorem, and the concept of normalization 
invoked is nonstandard. The properties of the balance function are consequently nonintuitive. 

• Cooling and diffusion time: According to arguments in [1] diffusion and cooling by isen- 
tropic expansion drive the charge-dependent correlation length in opposite directions. If the 
correlation length decreases it is therefore not clear what to conclude. Diffusion increases 
monotonically with the duration of hadronic rescattering. In [1] cooling by isentropic expan- 
sion occurs in the prehadronic phase. Cooling-induced narrowing would thus also signal a 
long-lived QGP. But this picture is at odds with a RHIC scenario of chemical decoupling at 
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~ 170 MeV followed by hadronic isentropic cooling and kinetic decoupling at ~ 120 MeV. 
Width reduction by hadronic cooling should signal an extended hadronic interval, not a long- 
lived QGP. 

• Balance function width and correlation length: The balance function was formulated by 
analogy with previous measures of correlation lengths in elementary collisions. It is thus 
implied that the balance function width is monotonically related to the charge-dependent 
correlation length on rapidity, if not equivalent to it. This implication is not valid. 



3 Charge correlations in elementary collisions 

A series of correlation measurements on rapidity in e^-e^, p-p and p-p collisions was critical to 
formulating the phenomenological representation of nonperturbative QCD contained in the Lund 
model, and for verifying that quarks carry charge and flavor quantum numbers. The experimental 
techniques developed in those analyses were invoked to motivate the balance function. It is useful 
to examine them in detail with reference to applications in A- A collisions. 



3.1 Techniques 

Correlations of conserved measures (charge, baryon number and strangeness) on azimuth and axial 
rapidity (or jet thrust axis in e~^-e~ collisions) were studied extensively for elementary collisions 
[4]. Several closely-related quantities were defined to measure two-particle rapidity correlations: 
associated charge density (PEP4), rapidity differences with tagged intervals (OPAL) and charge 
compensation probabilities (TASSO), all based on the conditional-probability concept. A schematic 
of a typical analysis is shown in Fig. 1. To accommodate limited statistics the two-particle or joint 
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Figure 1: Typical correlation anal- 
ysis of elementary collisions. uq 
is a kinematic limit. Gray bands 
represent long-range and short-range 
measure conservation on rapidity dif- 
ference t/A- Projections on one 
marginal axis are shown for condi- 
tions (A,B) on the other, with a car- 
toon of the system space-time evolu- 
tion. 



distribution on rapidity is typically divided into two regions by conditions placed on the rapidity 
of one (test or tagged) particle. The joint distribution for particles of opposite 'charge' (Q, B, S) is 
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then projected onto the rapidity for a second particle. These techniques can be compared with the 
structure impUed by Bayes' theorem as illustrated in Fig. 2 (left panel). The conditional projections 
on rapidity typically contain peaked features reflecting correlation structure on rapidity difference. 
The locations, widths, amplitudes and polarities of the peaks reveal long-range and short-range 
measure conservation among Q, B and S on a 'string' or ID QCD color field. 

3.2 Results 

The correlation structure typically has two peaked components, a short-range component (SRC) 
due to local measure conservation during string fragmentation, and a long-range component (LRC) 
representing the conserved quantum numbers of the primary quark or nucleon pair, still localized 
near the initial or limiting rapidities (±yo) in elementary collisions. The LRC reveals how the 
initial-state (q-q, p-p) quantum numbers (B, Q or S) propagate from the kinematic limit inward 
toward midrapidity (part of equilibration). The SRC reveals in part how intermediate-state quan- 
tum numbers created locally in pairs propagate outward from the formation region (also part of 
equilibration) . 

Correlations in elementary collisions revealed that a ID string has a natural correlation length 
for hadron pair correlations (SRC) ~ 1 fm {i.e., one hadron diameter) [5]. The fragmentation picture 

is consistent with hadrons initially forming a contiguous arrangement (a dense ID fluid). Observed 
short-range charge (baryon number, strangeness) correlations reflect hadron 'charge ordering' on a 
statistical basis: in the ID system it is more probable to find alternating charge signs, a form of 
local measure conservation. 

3.3 Implications for the balance function and A-A collisions 

Comparison of p-p collisions with A-A ccntrality dependence is critical to many aspects of heavy ion 
collisions, including charge-dependent correlations. According to [1] 'hadronization is nearly instan- 
taneous' in p-p collisions. Whether this means the interval from projectile contact to hadronization, 
the hadronization process itself or the time from hadronization to kinetic decoupling is not clear. 
The last is nominally accessible to the balance function. 

Assuming the last interpretation the SRC rapidity correlation length in p-p collisions should be 
determined by the thermal velocity distribution, the Bj or ken-expansion flow velocity and the char- 
acteristic length of ID spatial hadron correlations at hadronization; the balance function width in 
p-p should acquire a minimum value which A-A collisions with 'late' hadronization would approach 
asymptotically. In [1] is the statement '...compared to p-p collisions one expects the peak in the 
balance function in nucleus-nucleus collisions to be narrower....' which appears to be inconsistent. 

In [1] it is assumed that in A-A collisions the hadron pair-separation length scale at formation 
is zero, inconsistent with the fragmentation picture developed from elementary collisions. Any 
nonzero width of rapidity correlations should then be attributed only to diffusion and thermal 
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velocities. Hadron diffusion in elementary collisions is said to be negligible, yet the charge-dependent 
correlation length is expected to be larger, again seemingly inconsistent. 

Systematics of hadron production in elementary collisions provide some initial guidance for A-A 
collisions, but the presence of a dense medium and the dimensionality difference of the hadronization 
process may have dramatic impact. The string system with its axial symmetry yields a compact 
formulation of ID charge-correlation systematics. A similar description for hadronization in 2-3 
dimensions in A-A would require a more general analysis of charge correlations in 3D momentum 
space. Restricting to a single momentum component in A-A analysis by analogy with elementary 
collisions could be misleading. 

Invoking the balance function argument in its literal details we might predict an increase in 
balance-function width with increasing A-A centrality, from a p-p minimum to a maximum for mid- 
peripheral A-A due to increased diffusion time. An extended QGP duration in the more central 
events would then result in a reduction of the width toward the p-p limiting value. 

4 Balance Function Definitions 

The balance function is nominally defined as a linear combination of conditional probabilities: the 
probability that a particle of charge a found within a condition interval on rapidity is accompanied 
by a particle of charge b separated from a by a specified rapidity difference. I present three equivalent 
definitions, the first as presented in [1, 2], the second based on binning the single-particle rapidity 
distribution, the third in terms of properly-defined conditional probabilities. These definitions form 
the basis for algebraic study of the balance function. 

4.1 Theory definition 

This definition was formulated by analogy with analysis of elementary collisions in terms of joint 
and marginal rapidity distributions. Using the notation in [2], Ai,A2 in Eq. (1) are defined as 
conditions on some combination of the rapidity variables of the joint distribution 



where Ai : A2 is read pairs with condition Ai given that condition A2 is satisfied. Each A^fe(A2) 
or Nabi^i, ^2) is an ensemble-averaged event-wise number of particles of charge b or pairs of 
charge combination ab respectively, satisfying the rapidity combination (Ai,A2). The conditions 
are taken to be a rapidity acceptance and a rapidity difference respectively: A2 — > Ay, Ai k6y. 
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In particular, domains Ai are diagonal strip integrals on the joint rapidity distribution that differ 
markedly from the usual conditional probability definition. 

4.2 Definition based on uniform binning 

A binned rapidity axis y, with acceptance Ay, bin size 5y and bin number M(Ay, supports 
a particle distribution. Bin i has integrated multiplicity nai{5y) of 'charge' species a. The total 
multiplicity of charge a in the acceptance is J^^^f"^'^^^ ^aiiSy) = Na{Ay). The bin counts represent 
integrals of the form 

fiaiiSy) ~ / Pinaiy)dy (2) 
Jyi-Sy/2 

where number density pinaiu) represents a single-particle inclusive parent distribution estimated 
by the histogram of ensemble averages {fiaiiSy)} and from which each event provides a sample 
{nai{Sy)}. A similar description holds for pair counts from the two-particle distributions. The 
balance function histogram {B^} is then defined as 

+ J M{Ay,5y)-k 

Bk{Ay,Sy) = 1/2 ^ -a6 ^ 5Z riai ■ {nb(^i+k) - SabSko) A; G [0, M - 1] (3) 

a,b=- 1^1=1 ''^ai j=i 

+ ^ M{Ay,5y) 

= 1/2 ^ -aft _ Uai ■ {nb(i+k) - SabSko) ke[l-M,0] 

a,b=- 1^1=1 '^ai i=i-k 

where the delta functions represent omission of self-pairs. This definition is equivalent to those in 
[1, 2] but makes explicit the geometry of the strip integrals, and facilitates subsequent comparison 
with the autocorrelation distribution. 



4.3 Definition based on conditional probabilities 

Basic two-point and single-point probabilities can be defined in terms of joint and marginal bin 
multiplicities as Pij{ab) = nai{nbj - SabSij) / Na{Nb - Sab) and pi{a) = UaiiNb - dab)/Na{Ni, - Sab)- 
Bayes' theorem for ensemble means is then 



p,.(a6. Ay, Sy) . = "-^"'^ " (4) 

Pi{a) UaiiNb -Sab) 

This is the conditional probability that a particle of charge b will appear in bin j given that bin i is 
occupied by a particle of charge a as determined by the joint distribution Uaiijibj — SabSij), as shown 
in the left panel of Fig. 2, of which distribution nai{Nb — Sgj,) is a marginal projection. Elementary 
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collisions analyses adopted this conventional conditional-probability technique, projecting onto one 
of the marginal axes with a condition imposed on the other (typically dividing the accessible rapidity 
into two bins) as shown in Fig. 1. The result provides simple interpretation of projected peak widths 
in terms of correlation lengths. The cost of the method is a skewed projection of the underlying 
autocorrelation distribution onto the marginal axis, a nonoptimal compression technique. The 
technique is however commensurate with the low statistics available in early analyses. For A-A 
collisions a more precise treatment is justified. 




\ 


bj(S>')-5ab5,j) 
\ 


1 1 1 1 1 1 
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Figure 2: Standard con- 
ditional probability defini- 
tion (left panel) according to 
Bayes' theorem and definition 
associated with balance func- 
tion (right panel). 



The conditional geometry of the balance function is modified as shown in Fig. 2 (right 
p_anel), with the equivalences nai{ 5y) Na{Ay), naiimj - 5abSij) Nk{ab,5y) ^ iVafc(Ai, A2). 
Nk(ab) = J2i^ai ("-f)(i+fc) ~ SabSko) is a strip integral, a projection onto the difference variable j/a- 
The conditional probability according to the balance function treatment is then 



Pk:Ay{ab) = Nk{ab,Ay,6y)/Na{Ay) (iVb(Ay) - 6ab) 
With this definition the balance function takes the form 



(5) 



1 



Bk{Ay,6y) = -< E -abpk:Ay{ab) Na{Ay) {Nh{Ay) - 6ab)/NaiAy) 



, a, 6= 



= < Y.Pk:Ayiaa) - abpk:Ay{ab) NaiAy) Nh{Ay)/NaiAy) 



(6) 



with J2kPk:AyicLb, Ay,5y) = 1 for each ab combination by construction. 



5 The Autocorrelation and the Balance Function 



The autocorrelation, a special case of the cross correlation, is the projection of a two-point distri- 
bution onto its difference variables, a good compression strategy if correlations are approximately 
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space invariant (invariant on sum variable ys in Fig- 2) . It is a well-established correlation measure 
used in time-series analysis for many decades [6, 7]. The autocorrelation characterizes the internal 
correlations of a single distribution. The term 'cross-correlation' applies to mixed distributions, 
for instance of different charge types. For simplicity I refer to cross- and autocorrelations both as 
autocorrelations. A correlation length denoted by £,y typically refers to a characteristic length or 
width of the autocorrelation distribution (on the difference variable yA in Fig. 2). 



5.1 Autocorrelation definition 

In terms of 2D bin integrals of the joint distribution the autocorrelation for an aperiodic distribu- 
tion on a bounded domain [7] (also called an unbiased sample autocorrelation) applicable to e.g., 
distributions on y, rj is 

J M{Ay,Sy)-k 

Ak,obj{ab) = jj—j^ X] ■ {nh(i+k) - ^ab^ko) ke[0,M-l] (7) 

1=1 

J MiAy,5y) 

= M_^_J^ Yl ^"-i ' i'^Ki+k) - ^ab^ko) ke[l-M,0] 
i=l—k 

For simplicity in what follows I suppress expressions for negative values of k. The autocorrelation 
for periodic distributions on (f) is ylfc o5j(a6) = 1/M J2i^i ^^ai " (^^6(1+*;) ~ ^ab^ko)^ and is symmetric 
about the midpoints \k\ = {M — l)/2 for real nj. The autocorrelation (omitting self pairs) is 
thus composed of factorial moments and cross moments at the binning scale 5y, with the weighted 
integral 

M-l 

X {M - k)Ak,ob3{o.b) = Na{Ay) (NkiAy) - 5ab) (8) 
k=l-M 

It follows that 



Na{Ay) {NbiAy) - 6ab) Pk:Ay{ab) = (M - k)Ak,obj{ab) (9) 
5.2 Autocorrelation reference and net autocorrelation 

The differential or net autocorrelation depends on a reference definition. For this study I assume a 
factorization reference constructed from the inclusive single-particle distribution 

-j^ M-k 

Ak,ref{ab) = Jf—j^ Yl riaiUbii+k) (10) 
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with integral J2kLi-M(-^ ~ ^)^k,refi(^b) = Na{Ay) Nf,{Ay) and chargc-indcpcndcnt (CI) and 
charge-dependent (CD) linear combinations Ak^ref{CI) = Z)a6 ^k,refio-b) = J2iLi^ nifii+k ~ 
and Ak^ref{CD) = J2ab ah Ak^ref{ab) = Ei^i'' ^QinieQ(i+k)ni+k » 0, where bin-wise 

charge asymmetry eqi = Qi/fii = Ea=- «"ai/ Ea=- ^ai- 

The net autocorrelation is defined as AAk{ab) = Ak^obj{o,b) — Ak^ref{o,b) with CI and CD linear 
combinations. I decompose the CD net autocorrelation into its integral (the net-charge variance at 
scale 6y), and the zero-integral CD net autocorrelation AAk{CD) 



AAk{ab) = AAk{ab) + Ak,ref{ab) (11) 

^ Ak,refiCI) 



AAkiCD) ~ AAkiCD) + 

« AAkiCD) + a'^idy)-n/M 



iV(Ay) 



N{Ay) 



The zero-integral CD net autocorrelation AAk(CD) records charge-dependent pair transport within 
events, averaged over the ensemble. This decomposition clarifies the role of net-charge fluctuations 
in the CD net autocorrelation and the balance function. This particular form of net-charge variance 
results from the assumption of a factorization reference, in which case the variance agrees with the 
elementary statistical definition. This is usually not the best reference choice for a precision analysis. 

5.3 Relating the balance function to the autocorrelation 

The balance function is a hybrid between a conditional probability and an autocorrelation. Its 
definition uses the language of conditional probabilities but incorporates a form of projection onto 
a difference variable in the manner of an autocorrelation. The two distributions can be related 
algebraically according to the above definitions (Eqs. (3) and (7) for balance function and autocor- 
relation in terms of binning and Eq. (9) for the form factor) 

B,(A,M . 1/2 ± .a^iJ^Lim^ (12) 

a,b=- " 
(M -k)Ak,ref{CI) / NAk,ref{CD) NAAk{CD) \ 
\ Ak,ref{CI) Ak,ref{CI) j 

/i NAAkjCD) ] 

where AAk{C D) / Ak^ref{C I) <C 1 is a correlation measure per final-state hadron pair. The addi- 
tional factor N gives the correlation measure per charged hadron. The overall sign convention for 
the balance function is opposite to isospin convention. 
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The unit-integral first factor is approximated by the triangle {l — \k\/M)/M\i the single-particle 
dN/dr] distribution is nearly uniform within the acceptance. In [1] a width reference is suggested: 
'In a globally equilibrated system... there would exist no correlation between the balancing charges, 
and the... width of the balance function would then correspond to the extent of single-particle 
emission in momentum space.' We find that the reference distribution is actually a normalized 
projection of the d'^N / dyidy2 pair distribution (as seen through the detector acceptance) onto the 
difference variable yA ~ k5y. The triangular shape of this projection is a source of difficulty in 
interpreting the balance function. The first two terms in the second factor represent net-charge 
fluctuations (integral of the CD net autocorrelation) in the acceptance, with canonical suppression. 
The zero-integral last term represents CD charge-pair transport, which may be characterized by a 
correlation length if the net autocorrelation is peaked near A; = 0. 



5.4 Balance function integral 

From the definition in terms of conditional probabilities in Eq. (6) we obtain 



^Sfc(Ay,<5y) 

k 



1 



1 + iV,(Ay)iV,(Ay) 
9 



a,b= 



(13) 



N{Ay) 



which also follows directly from Eq. (12). The accuracy of the approximation depends on charge 
asymmetry eg <^1, a condition well satisfied for QED charge and central A- A collisions near 
midrapidity. For baryon number, strangeness or noncentral collisions this condition may not be 
satisfied, as noted in [2]. The balance-function integral thus depends on the acceptance fraction of 



^ uncorrelated 




Ax/Q 1 



Figure 3: Relationship between balance- 
function integral and net-charge fluctuations 
plotted vs acceptance fraction of 17 = in. 
Sloping line shows the effect of canonical sup- 
pression without charge correlations. Curve 
suggests the effect of additional charge corre- 
lations. 



the total 47r charge distribution in the form of canonical suppression of net-charge fluctuations and 
nontrivial charge correlations, as shown in Fig. 3. For a Air detector cJq = and X]fc -^fc = 1- 

The appearance of an extraneous variance in a correlation measure is an example of measure 
bias, a cross coupling of one aspect of correlation analysis with another, producing a measure that 
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is difficult to interpret. Equivalently, the balance function in Eq. (12) can be seen as an integral 
equation for the charge-dependent correlations represented by the zero-integral CD net autocorre- 
lation. Charge correlations enter twice: once in the integral ctq/N and once in the detailed shape 
of the pair-transport distribution /S.Ak. The latter is by itself an example of an unbiased measure 
which reveals the fundamental correlation length and amplitude in an unambiguous manner. 



5.5 Balance function reference 

The balance function distribution should be compared to a reference representing both absence of 
charge-dependent correlations and conservation of net charge within the detector acceptance. The 
factorization reference used for this autocorrelation (and hence balance function) illustration was 
adopted to simplify illustration. An alternative reference can be generated directly from data by a 
mixing procedure. 

Event mixing in this case involves charge mixing and/or rapidity mixing. Charge mixing does 
not respect conservation of net charge in 47r, which is a constraint for real data. The net-charge 
variance then takes the grand-canonical limit fTg/iV ~ 1, driving the integral of the balance function 
to zero as seen in Eq. (13) and Fig. 3. Forming mixed events by selecting random particles from 
an event ensemble also mixes charges and produces a reference that doesn't respect conservation 
of net charge. Rapidity mixing randomizes the location of existing charges within the acceptance, 
which reduces /S.Ai^{CD) to zero in Eq. (12). The resulting reference retains the correct form factor 
and net-charge fluctuation term. The algebraic form representing rapidity mixing is 

Bk,ref{^y,5y) p,^^^(CI) 1 1 - | (14) 



which can be compared to a balance function derived from unmixed real data to reveal the nontrivial 
charge correlations represented by N AA^i.{CD)/Ak^ref{CI)■ 



5.6 Pair-ratio distributions 

One can also generate ratios of object and reference pair distributions for each charge combination 
to access the CD net autocorrelation directly [8]. Event-wise ratios of sibling-pair bin counts to 
mixed-pair bin counts are constructed, forming sibling pairs within single events and mixed pairs 
by matching each object event to its nearest neighbor in a space spanned by global event properties 
such as total multiplicity [9] and then averaging the event-wise bin ratios over an event class 



rij{ab) = {uai (ubj - SabSij)}sib / {nai {ribj - SabSij)}. 



(15) 
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The relative net autocorrelation is expressed in terms of this ratio distribution by 

1 



AAk{ab)/Ak,refiab) = 



M-k 



M(Ay,5y)-k 
i=l 



1 



(16) 



Because the sibling and mixed pair totals are matched on an event-wise basis the net autocorrelation 
AAk{ab) has zero integral by construction for each ab combination. This method accesses the 2D 
correlation structure of each charge combination, not only the CD or CI combination. The CD 
linear combination reveals the charge-dependent correlations which are the nominal object of the 
balance function. The pair-ratio method is acceptance independent (whereas ctq and B are not) 
and provides direct access to the basic correlation length and amplitude oq, as well as other 
details of CD pair transport. 



6 Balance Function Width Systematics 

The theoretical hypothesis which motivates the balance function is expressed in terms of parametric 
variation of a charge-dependent correlation length. The balance function is said to represent this 
correlation length through its width. I examine the relationship between balance function width and 
correlation length with a simple model that corresponds closely to correlation structure observed 
in A- A collisions. 

6.1 Correlation model 

The model is formulated in terms of ratio AAj.{CD)/Ak^ref{CI) derived from a nominally uniform 
distribution of charges on momentum variable y (which could be rapidity y, pseudorapidity rj or 
azimuth (p). Given difference variable j/a = yi — 2/2 I define normalized difference variable x = 
yA/Ay G [—1, 1] and normalized correlation length xq = ^y/Ay G [0, oo), both relative to detector 
acceptance Ay. This continuum model relates to the discrete case by k/M kSy/MSy ^ x. I 
adopt a simple exponential form for the relative net autocorrelation consistent with observations: 
—AA{x)/A(x) = ao {e~'^l/^° — Aq}. The integral condition 

j\l-\x\)AA{x)/A{x)dx = (17) 

is imposed to determine the offset Aq, from which we obtain 

Ao = 2x0 {1 - ^0 (1 - e-V^o)} (18) 

The net-autocorrelation model is shown in Fig. 4 (right panel) for a selection of eight relative 
correlation lengths xq plotted vs rapidity difference assuming Ay = 2, with amplitude ao held fixed 
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Figure 4: Autocorrelation distribu- 
tions (right panel) for a selection of 
eight correlation lengths with ampli- 
tude held fixed [8] and correspond- 
ing distribution values (left panel) at 
X = 1 (two correlation lengths are 
outside the plot range). 



and offset Aq varying according to Eq. (18) to maintain the zero- integral condition. Distribution 
values at yA = 2 or x = 1 are shown in the left panel to illustrate the systematic variation of Aq. 
Defining Aq = N qq and Qq = Uq/N, the balance function for this idealization is then 



(l-|x|){l-Qo + ^o (e-N/^°-Ao)} 



with the condition J^^ B{6y,x) dx = 1 — Qo^Ay). The basic correlation parameters and are 
best obtained directly from the CD net autocorrelation, from which we can derive other related 
quantities, in particular the balance function. 

6.2 Width systematics 

We now investigate the parametric dependence of the balance function width on the model pa- 
rameters - the amplitude and correlation length. I define the mean-square width of the balance 
function as 



1 

— / x^B{x)dx (20) 

1- Qo J -I 



1 ^ 2Ao 



6 1 - Q, 



{2x1 - 6^) (1 - e-'/^o) + {6x1 + 4) - 3^ Ao 



with ay = Ay ■ ax being the corresponding rms width on rapidity difference. The rms width for 
the case of no correlations (^o = 0) is ^1/6 ~ 0.4 times Ay. For the STAR detector the latter is 
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Figure 5: Balance functions (left 
panel) corresponding to eight auto- 
correlations in Fig. 4 and distri- 
butions of rms widths (right panel) 
for four amplitudes (^o) in Eq. (20). 
The dots correspond to balance func- 
tions in the left panel (two outside 
the plot range). 



about 2.1, and the uncorrelated or rapidity-mixed reference rms is observed to be about 0.83 as 
shown in Fig. 5 (right panel). 

The balance functions in the left panel exhibit the expected behavior: charge correlation 
structure deviating from a triangular reference distribution whose area is determined by net-charge 
fluctuations within the acceptance. The surprising result is in the right panel. The horizontal line 
at 0.83 approximates the reference rms width corresponding to no charge correlations. The balance 
function rms width otherwise has a minimum for ^y/Ay « 1/3, which coincides with typical values 
for heavy ion collisions and the STAR detector (lower horizontal bar) [8]. For relative correlation 
lengths much below this value the width of the balance function would actually anticorrelate with 
the charge-dependent correlation length. For heavy ion collisions the balance function rms width 
is thus dominated by the correlation amplitude, as illustrated by the four curves in the right panel. 
It is essentially insensitive to the charge-dependent correlation length. 

7 Conclusions 

The balance function, proposed as a method to probe QGP existence and duration by means 
of charge-dependent correlations, has stimulated re-examination of local measure conservation in 
heavy ion collisions and its analysis methods. I present here several critical issues concerning the 
motivating chain of argument and interpretation of the balance function as a correlation measure. 

If the charge-dependent net autocorrelation has a simple peaked structure near zero momen- 
tum difference ('short-range' correlations, as observed in both elementary and A-A collisions) the 
basic charge-dependent correlation parameters are the peak amplitude and its rms width, the latter 
conventionally called the correlation length. Detailed aspects of the distribution are also impor- 
tant, but for the balance function hypothesis these are the relevant quantities. The purpose of a 
correlation analysis applicable to the physics of the balance function should be inference of this 
correlation length and amplitude. 

The balance function apparently confuses amplitude and correlation length in its width defini- 
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tion. For heavy ion collisions viewed by a large-acceptance detector (STAR at RHIC) the balance 
function width is ironically insensitive to the correlation length, is dominated instead by the cor- 
relation amplitude. Theoretical hypotheses addressing the charge-dependent correlation length in 
terms of the balance function width are consequently misdirected. 

The supporting theoretical chain of argument concerning initial charge correlations at hadroniza- 
tion, exclusive influence of diffusion and thermal velocity on final-state correlation lengths and 
relation of A- A correlations to those in elementary collisions seems iconsistent. The structure of 
hadronization in A-A collisions is an open question; extrapolations from ID string phenomenol- 
ogy are unwarranted without careful re-examination. The geometry of the prehadronic medium at 
hadronization may dominate the structure of A-A charge correlations. 
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